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Abstract 

We explore a connection of the forced Burgers equation with the Schrddinger 
(diffusive) interpolating dynamics in the presence of deterministic external 
forces. This entails an exploration of the consistency conditions that al- 
low to interpret dispersion of passive contaminants in the Burgers flow as 
a Markovian diffusion process. In general, the usage of a continuity equation 
dtp = —V(vp), where v = v(x,t) stands for the Burgers field and p is the 
density of transported matter, is at variance with the explicit diffusion sce- 
nario. Under these circumstances, we give a complete characterisation of the 
diffusive matter transport that is governed by Burgers velocity fields. The 
result extends both to the approximate description of the transport driven by 
an incompressible fluid and to motions in an infinitely compressible medium. 
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The Burgers equation with (typically without, |T], |2|) the forcing term F(x,t): 



d t v + (vV)v = vAv + F(x,t) (1) 

and especially its statistically relevant curlv = solutions (in below we shall not 
use random initial data), recently have acquired a considerable popularity in various 
physical contexts, P]-[]r9|. 

Although Burgers velocity fields can be analysed on their own, frequently one 
needs a supplementary insight into the matter transport dynamics that is consistent 
with the chosen (Burgers) velocity field evolution. Then, the passive scalar (tracer 



or contaminant) advection-in-a-flow problem, [14, 11, 16 naturally appears through 
the parabolic dynamics: 

d t T + (tJV)T = uAT (2) 

While looking for the stochastic implementation of the microscopic (molecular) dy- 
namics (2), |2T| , |Tl~| , |16|j , it is assumed that the "diffusing scalar" (contaminant in 



the lore of early statistical turbulence models) obeys an Ito equation: 

dX(t) = v(x, t)dt + V2udW{t) (3) 

X(0) = f _> X(t) = x 

where the given forced Burgers velocity field is perturbed by the noise term rep- 
resenting a molecular diffusion. In the Ito representation of diffusion-type random 
variable X{t) one explicitly refers to the Wiener process \^2vW(t), instead of the 
usually adopted formal white noise integral / 'f/(s)(is, coming from the Langevin- 
type version of (3). 

Under these premises, we cannot view equations (l)-(3) as completely indepen- 
dent (disjoint) problems: the velocity field v cannot be arbitrarily inferred from (1) 
or any other velocity-defining equation without verifying the consistency conditions, 
which would allow to associate (2) and (3) with a well defined random dynamics 
(stochastic process), and Markovian diffusion in particular, |2~3fl . 

In connection with the usage of Burgers velocity fields (with or without external 
forcing) which in (3) clearly are intended to replace the standard forward drift of the 
would-be-involved Markov diffusion process, we have not found in the literature any 
attempt to resolve apparent contradictions arising if (2) and/or (3) are defined by 
means of (1). Also, an issue of the necessary correlation (cf. |Tj|, Chap. 7. 3, devoted 
to the turbulent transport and the related dispersion of contaminants) between the 
probabilistic Fokker-Planck dynamics of the diffusing tracer, and this of the passive 
tracer (contaminant) concentration (2), has been left aside in the literature. 
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Moreover, rather obvious hesitation could have been observed in attempts to 
establish the most appropriate matter transport rule, if (l)-(3) are adopted. De- 
pending on the particular phenomeno logical departure point, one either adopts the 
standard continuity equation, || ||, that is certainly valid to a high degree of ac- 
curacy in the low viscosity limit v J, of (l)-(3), but incorrect on mathematical 
grounds if there is a diffusion involved and simultaneously a solution of (1) stands 
for the respective current velocity of the flow: d t p(x,t) = — V[v(x, t)p(x, t)] . Al- 
ternatively, following the white noise calculus tradition telling that the stochastic 
integral X(t) = J * v(X(s), s)ds + fQff(s)ds necessarily implies the Fokker-Planck 
equation, one adopts, pif : d t p(x,t) = vAp(x,t) — V[v(x, t)p(x, £)] which is clearly 
problematic in view of the classic Mc Kean's discussion of the propagation of chaos 
for the Burgers equation, E|, I2R I2H] and the derivation of the stochastic " Burgers 



process" in this context: "the fun begins in trying to describe this Burgers motion 
as the path of a tagged molecule in an infinite bath of like molecules" , [B3J] . 

To put things on the solid ground, let us consider a Markovian diffusion process, 
which is characterised by the transition probability density (generally inhomoge- 
neous in space and time law of random displacements) p(y, s, x, t) , < s < t < T, 
and the probability density p(x,t) of its random variable X(t) , < t < T. The 
process is completely determined by these data. For clarity of discussion, we do not 
impose any spatial boundary restrictions, nor fix any concrete limiting value of T 
which, in principle, can be moved to infinity. 

The conditions valid for any e > 0: 

(a) there holds Um tls -^ !\$-$\ >e P(y\ s > t ) ( ^ x = °> 

(b) there exists a (forward) drift b(x,s) = Um t i s ^ ij^_^i< e (y — x)p(x, s,y,t)d 3 y, 

(c) there exists a diffusion function (in our case it is simply a diffusion coefficient v) 
a(x, s) = Um tis ^- s J w ^\< £ (y - x) 2 p(x, s,y,t)d 3 y, 

are conventionally interpreted to define a diffusion process, PB] , Under suitable 
restrictions (boundedness of involved functions, their continuous differentiability) 
the function: 

g(x, s) = J p(x, s, y, T)g(y, T)d 3 y (4) 

satisfies the backward diffusion equation (notice that the minus sign appears, in 
comparison with (2)) 

-d s g(x, s) = vAg(x, s) + [b(x, s)V]g(x, s) . (5) 

Let us point out that the validity of (5) is known to be a necessary condition for the 
existence of a Markov diffusion process, whose probability density p(x, t) is to obey 
the Fokker-Planck equation (the forward drift b(x, t) replaces the previously utilized 
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Burgers field v(x,t))): 

8 t p(x, t) = uAp(x, t) - V[6(x, t)p(x, t)] (6) 

The case of particular interest in the nonequilibrium statistical physics literature 
appears when p(y, s, x, t) is a fundamental solution of (5) with respect to variables 
y, s, ^2], E3[], see however |28| for an alternative situation. Then, the transition 
probability density satisfies also the second Kolmogorov (e.g. the Fokker- Planck) 
equation in the remaining x, t pair of variables. Let us emphasize that these two 
equations form an adjoint pair , referring to the slightly counterintuitive for physi- 
cists, although transparent for mathematicians, [31], [H| |34j| , issue of time 
reversal of diffusions. 

After adjusting (3) to the present context, X(t) = J b(X(s), s) ds + y/2uW(t) 
we can utilize standard rules of the Ito stochastic calculus, [^, [32], |3^, to realise 
that for any smooth function f(x,t) of the random variable X(t) the conditional 
expectation value: 

lim mQ ^ t [j p(x, t,y,t+ At)f(y, t + At)d 3 y - f(x, t)} = (D+f)(X(t),t) = (7) 

= (d t + (W) + vA)f(x,t) 
X(t) = x 

determines the forward drift b(x, t) (if we set components of X instead of /) and 
allows to introduce the local field of forward accelerations associated with the dif- 



fusion process, which we constrain by demanding (see e.g. Refs. [[35], [33], [34J for 
prototypes of such dynamical constraints): 

(D 2 + )X)(t) = (D + b){X(t),t) = (d t b+ (bV)b+ uAb)(X(t),t) = F(X{t),t) (8) 

where, at the moment arbitrary, function F(x, t) may be interpreted as the external 
forcing applied to the diffusing system, [p5[] . In particular, if we assume that drifts 
remain gradient fields, curlb = 0, under the forcing, then those that are allowed by 
the prescribed choice of F(x, t) must fulfill the compatibility condition (notice the 
conspicuous absence of the standard Newtonian minus sign in this analogue of the 
second Newton law) 

F(x , t) = VQ{x , t) (9) 
Sl(2,t) = 2v[dt* + ~(|; + V&)] 
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establishes the Girsanov-type martingale connection of the forward drift b(x, t) = 
2z/V$(x, £) with the (Feynman-Kac, cf. [29, 28[) potential Q(x, t) of the chosen 
external force field. 

One of distinctive features of Markovian diffusion processes with the positive 
density p(x, t) is that the notion of the backward transition probability density 
p*(y, s, x,t) can be consistently introduced on each finite time interval, say < 
s < t < T: 

p(x, t)p*(y, s, x, t) = p(y, s, x , t)p(y, s) (10) 

so that / p(y,s)p(y,s,x,t)d 3 y = p(x,t) and p(y,s) = J p*(y, s, x,t)p(x,t)d 3 x. This 
allows to define the backward derivative of the process in the conditional mean (cf. 



35, 29, 36, p% for a discussion of these concepts in case of the most traditional 



Brownian motion and Smoluchowski-type diffusion processes) 

Um Atl0 -^[x - J p*{y,t - At,x,t)yd 3 y] = {D^X){t) = t{X{t), t) (11) 

(D_f)(X(t),t) = (d t + (6*V) - uA)f(X(t),t) 

Accordingly, the backward version of the dynamical constraint imposed on the ac- 
celeration field reads 

(D*_X)(t) = (DlX)(t) = F(X(t),t) (12) 

where under the gradient-drift field assumption, curl b* = we have explicitly ful- 
filled the forced Burgers equation (cf. (1)): 

d t t + (KV)t - vAb* = F (13) 



where, [35], [33], p9l , in view of 6* = b — 2uVlnp, we deal with F(x,t) previously 
introduced in (9). A notable consequence of the involved backward Ito calculus is 
that the Fokker-Planck equation (6) can be transformed to an equivalent form of: 

d t p{x, t) = -uAp{x, t) - V[K(x, t)p(x, t)] (14) 

with the very same initial (Cauchy) data po(x) = p(x,0) as in (6). 

At this point let us recall that equations (5) and (6) form a natural adjoint 
pair of equations that determine the Markovian diffusion process in the chosen time 
interval [0, T]. Clearly, an adjoint of (14), reads: 

d s f(x, s) = uAf(x, s) - [l(x, s)V]f(x, s) (15) 
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where: 

f(x, s) = I p*(y, 0, x, s)f(y, 0)d s y , (16) 



to be compared with (4), (5) and the previously mentioned passive scalar dynamics 
(2), see e.g. also [^TJ. Here, manifestly, the time evolution of the backward drift is 
governed by the Burgers equation, and the diffusion equation (15) is correlated (via 
the definition (10)) with the probability density evolution rule (14). 

This pair only can be consistently utilized if the diffusion proces is to be driven 
by forced (or unforced) Burgers velocity fields. 

Let us point out that the study of diffusion in the Burgers flow may begin from 
first solving the Burgers equation (12) for a chosen external force field, next specify- 
ing the probability density evolution (14), eventually ending with the corresponding 
"passive contaminant" concentration dynamics (15), (16). All that remains in per- 
fect agreement with the heuristic discussion of the concentration dynamics given in 
Ref. [IE||, Chap. 7.3. where the "backward dispersion" problem with "time running 



backwards" was found necessary to predict the concentration. 

Let us notice that the familiar logarithmic Hopf-Cole transformation, [^|, |3"5| , 
of the Burgers equation into the generalised diffusion equation (yielding explicit 
solutions in the unforced case) has received a generalisation in the framework of the 
so called Schrodinger boundary-data (interpolation) problem, 0, ^9], [37|, |38| , 



see also |j40| , (0J. In particular, in its recent reformulation, |29|, P8| , the problem 
of defining a suitable Markovian diffusion process was reduced to investigating the 
adjoint pairs of parabolic partial differential equations, like e.g. (5), (6) or (14), (15). 
In case of gradient drift fields this amounts to checking (this imposes limitations on 
the admissible force field potential) whether the Feynman-Kac kernel 

k(y,s,x,t) = J exp[- J c(w(T),r)dr]d^jj(w) (17) 

is positive and continuous in the open space-time area of interest, and whether it 
gives rise to positive solutions of the adjoint pair of generalised heat equations: 

dtu(x, t) = isAu(x, t) — c(x, t)u(x, t) (18) 

d t v(x, t) = —uAv(x, t) + c(x, t)v (x, t) 



(x,t) 
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where c(x, t) = ^Vt{x, t) follows from the previous formulas. In the above, d/J,% 
is the conditional Wiener measure over sample paths of the standard Brownian 
motion. 

Solutions of (18), upon suitable normalisation give rise to the Markovian diffusion 
process with the factorised probability density p(x,t) = u(x,t)v(x,t) wich interpo- 
lates between the boundary density data p(x, 0) and p(x,T), with the forward and 
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backward drifts of the process defined as follows: 



b(x,t) = 2v^£ (19) 
v(x, t) 

7* /- \ V«(f, t) 

u(x, t) 

in the prescribed time interval [0, T\. The transition probability density of this 
process reads: 

p(y, s, x, t) = k(y, s, x, t) V X J \ , (20) 

Here, neither k, (17), nor p, (20) need to be the fundamental solutions of appropriate 
parabolic equations, see e.g. ref. |28| where an issue of differentiability is analyzed. 

The corresponding (since p(x,t) is given) transition probability density, (10), of 
the backward process has the form: 

u(y^ s) 

P* (y, s, x, t) = k(y, s, x, t) ' . (21) 

u{x, t) 

Obviously, |28| , [33|| , in the time interval < s < t < T there holds: 
u(x, t) = J u (y)k(y, s, x, t)d 3 y and v(y, s) = J k(y, s, x, T)v T {x)d?x. 

By defining = logu, we immediately recover the traditional form of the 
Hopf-Cole transformation for Burgers velocity fields: = — 2iA7$*. In the special 
case of the standard free Brownian motion, there holds b(x,t) = while b*(x,t) = 
—2vVlog p(x, t). 
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